Subgap tunnelling through channels of polarons and bipolarons in chain conductors. 
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(Dated:) 

We suggest a theory of internal coherent tunnelling in the pseudogap region where the applied 
voltage is below the free electron gap. We consider quasi ID systems where the gap is originated 
by a lattice dimerization (Peierls or SSH effect) like in polyacethylene, as well as low symmetry ID 
semiconductors. Results may be applied to several types of conjugated polymers, to semiconducting 
nanotubes and to quantum wires of semiconductors. The approach may be generalized to tunnelling 

» ' in strongly correlated systems showing the pseudogap effect, like the family of High Tc materials in 

, the undoped limit. We demonstrate the evolution of tunnelling current-voltage characteristics from 

smearing the free electron gap down to threshold for tunnelling of polarons and further down to the 

I ' region of bi-electronic tunnelling via bipolarons or kink pairs. The interchain tunnelling is described 

in a parallel comparison with the on chain optical absorption, also within the subgap region. 

, PACS numbers: 05.60.Gg 71.10.Pm 71.27. +a 71.38.-k 71.38.Mx 73.40.Gk 78.67.-n 79.60.-i 

£j ' I. INTRODUCTION. 

r/3 ; 

The interchain, interplane transport of electrons in low dimensional (quasi ID, 2D) materials attracts much 
attention^ in view of striking differences between longitudinal and transverse transport mechanisms revealing a gen- 
eral problematics of strongly correlated electronic systems. Beyond the low field (linear) conduction, the tunnelling 
current- voltage J-U characteristics J(U), a = dJ/dll are of particular importance. The interest has been renewed 
thanks to recently developed^ design of intrinsic tunnelling devices where electronic transitions between weakly coupled 
O \ chains or planes take place in the bulk of the unperturbed material. 

i ^ i ■ The first feature one expects to see at any tunnelling experiment in gapful conductors is the regime of free electrons 
when the current onset corresponds to the voltage U = Eg of the gap in the spectrum of electrons. But contrarily to 
■ usual systems, like semi- or even superconductors, there is also a possibility for tunnelling within the subgap region 
Eg < U < Eg. It is related to the pseudogap (PG) phenomenon known for strongly correlated electrons in general, 
well pronounces in quasi ID systems and particularly in cases where the gap is opened by a spontaneous symmetry 
breaking (see^ and refs. therein). The PG is originated by a difference, sometimes qualitative, between three forms 
of electronic states: a) short living excitations which are close to free electrons, b) dressed stationary excitations of 
i— ( ■ the correlated systems, and c) added particles which modify the ground state itself. For our typical examples of 
\ electrons on a flexible lattice, the modification results in self-trapped states (b) like single particle v = 1 polarons 
with energies W\ < Ao below the single electron (a) activation energy Ao; then the new gap E g < Eg = 2Aq will be 
observed as a true threshold with the PG in between. There may be also contributions of two-particle v — 2 states (c) 
- bipolarons, which energy gain per electron is larger than for polarons Wi < 2W\. While the cases (a,b) are common 
for low symmetry and discrete symmetry cases, for (c) there is a further drastic effect of a spontaneous symmetry 
breaking like the case of the polyacethylene (CH) X or of some doubly commensurate CDWs. Now the bipolarons 
are decoupled into particles with a nontrivial topology, solitons or kinks, changing the sign of the order parameter of 
the dimerization. The situation is further intricate in systems with a continuous GS degeneracy like Incommensurate 
Charge Density Waves (ICDW) or Wigner crystals. Here even the self-trapping of a single electron is allowed to 
lead to topologically nontrivial states, the amplitude solitons ASs. In the same class we find a more common case of 
acoustic polarons in a ID semiconductor^^. 

Properties of systems with different types of the GS degeneracy, and required theoretical approaches, are quite 
different. Here we shall concentrate on systems with a discrete, precisely double, degeneracy which also include most 
basic elements of non-degenerate systems. Theoretically, the tunnelling in CDWs was studied in details for regimes of 
free electrons^ when the current onset corresponds to the voltage U = Eg of the gap in the spectrum of electrons. We 
shall consider the tunnelling in the PG regime. We shall follow the method^ developed for studies of single particle 
spectral density I(p,lu) in applications to PES and ARPES intensities. We refer to this publication for details in 
techniques and literature. 

A word of notations. In the following we shall invoke many quantities with the dimension of energy (or frequency, 
since we shall keep h = 1) which will be classified according to different characters (with indices). U > and 
fl > will be the external voltage difference for tunnelling and the external frequency for PES or optics. E > 
will always stay for electronic eigenvalue in a given potential, negative values will be addressed explicitly as —E. 
V U (E) will be branches of a total energy (of deformations together with electronic energies) supporting eigenstates 
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±E which may be filled with occupation numbers v = 0, 1, 2. W v will be total energies of stationary states, that is 
W v — min s V V (E) = V V {E V ). uj <C E will be the frequency of a collective mode (phonons specifically to CDWs) which 
interaction with electrons is responsible for their self-trapping. The collective deformation A(x) will also be measured 
as the potential energy experienced by electrons. We shall keep the electron charge e = 1 hence potentials will be 
measured as energies and the interchain current J will have the dimension number- of-particles /unit-time /unit-length. 
The indices j = a,b will number coupled chains; indices i will number moments r, of time for virtual processes. 

II. SPECTROSCOPIES OF THE PSEUDOGAP. 

A possibility of tunnelling or of other excitations within the gap in spectra of free electrons E < Aq is related 
to a more general phenomenon of the pseudogap PG. For electrons, the PG signifies the remnants of the spectral 
density P), or the integrated one J(Q) = J I(Q, P)dP/2n, at Wi < il < Ao where Wi is the absolute boundary 
of the spectrum. W\ is the energy of a fully dressed state of one electron interacting with other degrees of freedom. 
(There may be totally external modes like deformations or polarizations for usual polarons, external modes essentially 
modified by the bath of electrons like in CDWs, internal collective modes of electronic system itself like in SDWs.) 
Most commonly, the self-trapped state of one electron is known as the "polaron " while more complex objects, solitons, 
can appear for systems with continuously degenerate GSs (see^ for a review). 

The functions and I(f2, P) are measured directly in PES and ARPES experiments (these abbreviations stay for 
the integrated Photo-Emission Spectroscopy and for the Angle (that is momentum) Resolved one) . As such they have 
been studied theoretically for the PG region by the present authors2*£*& and we refer to these publications for a more 
comprehensive discussion and for the literature review. The one electron spectra can be accessed also in traditional 
external tunnelling experiments: junctions or STM. For the last case, and practically for macroscopic point junctions, 
only the integrated J(Q) is measured. Elements of the full dependence J(f2,P) become necessary to describe strongly 
anisotropic materials (layered quasi-2D or chain quasi-lD ones) where the coherent tunnelling is realized in internal 
junctions of "mesa" type devices 2 . Here the tunnelling goes between adjacent layers within the single crystal of the 
same material, hence the momentum is preserved. In a simpler version, the internal subgap tunnelling takes place 
from free electrons of some metallic bands or pockets to polaronic states within gapful spectra which probably takes 
place in NbSey^. Otherwise it measures actually the joint spectral density for creation a particle-hole excitation at 
adjacent chains (the interchain exciton). In this respect it will be instructive to compare the coherent tunnelling and 
the subgap optical absorption OA (see a short excursion and references in£, III.E). A less expected version of the 
internal tunnelling is a possibility for bi-electronic transfers (tunnelling of bipolarons or of kink-antikink pairs) which 
usually is attributed only to superconductors. We shall see that these processes extend the PG further down to even 
lower voltages. 

In any case, the tunnelling current J(U) is given by the transition rate of electrons between two subsystems a, b kept 
at the potential difference U. For a weak coupling t±, the electron tunnelling from a to b is given by the convolution 
of spectral densities 

j~t\ i a (-n,p)i b (u -n,p)\A(P)\ 2 dsidP 

Jo 

if the momentum is conserved, or of their integrals J ~ t\ I a (—Q) Ib{U — Q)dfl for the incoherent tunnelling. 
(Everywhere we assume T = 0.) Recall that for free electrons with a spectrum E(P) we have 7(f2, P) ~ S(\Cl\ — E(P)) 
while I(Q) becomes the DOS N{Q), e.g. for D = 1 7(fi) - t^m 1 / 2 ^ - Ao)- 1 / 2 near the bottom of the free 

band where the electron effective mass is m. 

Consider briefly the case where one of reservoirs, say #6, is composed by free electrons with a known DOS Nb(Cl). 
One of applications of a tunnelling between the free spectrum and the PG may be the case of several families of 
conjugated polymers (polypyrolle, polythiophcnc) where origins of filled, ir or empty, ir* bands below and above 
the gap are essentially different. Then the polaronic effect, hence the PG, may be pronounced only for one type of 
particles: electrons or holes. The same concerns ID systems made of semiconducting wires where both effective masses 
and deformation potentials for electrons and holes are usually very different. Then, for the incoherent tunnelling, 
I a (U) gives directly either the tunnelling current J{U) ~ Ia(U) (if Nb has a sharp peak at the Fermi level, which is 
typical for using junctions with superconductors) or the tunnelling differential conductance a = dJ/dU ~ Ia(U) (if 
Nb rj const at the Fermi surface). 

Below we shall be mostly interested in systems with the charge conjugated symmetry (or qualitatively equivalent 
ones); the examples are carbon nanotubes, symmetric conjugated polymers like polyphcnylcncs, polyanilincs and 
polymers where the gap is formed (partly at least) by the spontaneous symmetry breaking: the polyacethylenes (— ). 
Numerical details will be presented for the last rich case. In all these cases the PG will exist near both rimes ±Ao of 
the free excitation gap E® = 2Aq. 
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1: The logarithmic plot for tunnelling or absorption intensities In / ~ — S in the pseudogap regime: between 2W\ and 



Recall now some known results for within the PGi It has the form / = Aexp(— S) where the action S = S(Q) 
is proportional to the big parameter of our adiabatic approximation: S ~ Ao/^o ^ 1- S(Q) is determined by an 
optimal fluctuation localized in space and time (an instanton) which supports the necessary split-off local level E. In 
principle, the prefactor A = A(f2) also depends on Q and may show power law dependencies near extremals 0, Wi. 
But within constraints of the adiabatic approximation 5£l too the dependence A(Q) is negligible in comparison with 
the one of 5(0). The characteristic value of A may be important for estimates of the overall magnitude of observable 
effects. Thus for the single particle integrated intensity A ~ (wWi) -1 ' 2 and Ap ~ (m/S.QUjWi)~ 1 / 2 for the momentum 
resolved intensity; here m is the effective electron mass m ~ Aq/v f . Appendix IbI contains derivation of the prefactor 
specifically for the tunnelling processes. 

In limiting cases we have 3 - 

1. Near the entry to the PG, just below the free edge Aq: 



A - W\ » A - il > : I = A exp 



enst 



Wo 



Ao 



3/2' 



2. Near the low end of the PG, just above the true spectral boundary W\: 



Ao - Wi > O - Wi > : I = Aexp 
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The total dependence and the values of numerical constants in the above limiting laws, can be determined 

approximately^ with the help of the instanton techniques simplified by the zero dimensional reduction (the Anzats of 
an effective particle which we shall recall and extend below). The resulting curve is plotted at the figure^ Moreover, 
the regime (I.) can be mapped exactly 3 " upon the problem of a particle in a quenched random uncorrelated potential 
which here is created by instantaneous quantum fluctuations of the media. The known exact solution^ provides the 
reference value of the coefficient in the exponent of Q, from which our approximate value differs only by 8% 3 . 

Recall for comparison the usual regime U > Eg of the allowed tunnelling which is dominated by free electronic 
states. The current of the coherent tunnelling between chains a, b is given as 



/ £<*M 



U - E a {p)) \K{p)\ 2 , A(p) = / dxVbpix)** (x) 



(3) 



with E a fi{p) f« ± (Ao +p 2 /2m e ) and ^j P (x) being the Bloch functions. It is instructive to compare the interchain 
tunnelling probability with the on chain interband optical absorption OA when the matrix element of density A(p) 
changes to the one of the momentum: t±A(p) Aoa(j>)- In both cases the e — h pair is created and the same spectral 
densities are involved. The difference is in matrix elements: the OA takes place between states of opposite parity 
while the tunnelling requires for the same parity. The on-chain OA between edges ±Ao of the free gap is known to 
be allowed since the parity of states near ±Aq is opposite, hence Aqa is finite and the OA intensity as a function of 
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FIG. 2: The logarithmic plot for bi-electron tunnelling or for optical absorption intensities In I ~ — S in the two particle 
pseudogap regime: between 2W 3 and 2Aq. Notice the ■\/Q — 2E S singularity near the two soliton threshold. 



frequency O rises as Ioa oc \j\JTT- ~ 2Aq. But for the same reason, the tunnelling matrix element between identical 
chains is prohibited at p — and the tunnelling will show only a weak edge onset J ~ — 2Arj. Nevertheless, in 
many cases of gaps opened due to spontaneous dimerization, the neighboring chains tend to order in antiphase. Now 
the shift by half a period permute states with E ^ then the parity of states near opposite rims ±Ao at neighboring 
chains is equal, the tunnelling becomes allowed and the usual singularity is restored: J oc l/\/U — 2Aq. 

Going down into the PG O < Ao, the above analysis applies to the on-chain optics but changes drastically for the 
interchain tunnelling. The tunnelling will be studied in details below, here we shall only mention in advance an effect 
of spatial incoherence of optimal quantum fluctuations at different chains which removes completely the constraints of 
orthogonality. The case of the on-chain OA can be analyzed briefly already here. The OA is given by the convolution 
of two fast decaying functions of the energy 

Ioa ~ J m u P)I(n - Oi, P) \^oa{P)\ 2 d^dP 

J exp [-5(Oi, P) ~ 5(0 - O x , P)] |A(Oi, O x - 0)| 2 dOi ~ 

A 2 (0/2, -0/2) (S"y 1/2 exp [-25(0/2)] = 

(S"(tt/2)r 1/2 / 2 (0/2); S" = 0. (4) 

Here we have used that for the convex function 5(0), as given by (11I2[I . the minimum of the expression 5(0i, P) + 
5(0 — Oi , P) lies at the middle Oi = U/2. At this point the electron levels E and E — O are placed symmetrically, 
wave functions have opposite parity, hence A(E, —E) 7^ is finite. This is the case of typical Peierls insulators. But 
for systems where the basis wave functions of valent and conductive bands have the same parity (the dipole OA is 
not allowed), Aoa(E, — E) = and we have to consider in the deviations from the symmetry condition. Now 
AoA(E a -, Eb) ~ (E a + Ei,) 2 and the saddle point integration in Q gives another factor of 1/5" which is small as 
~ lu /A . We arrive at the answer similar to 10} but with the small prefactor (S"')~ 3 ^ 2 . 

Until now we did not consider the dependencies on the momentum P. In the full range of O and P, the spectral 
function 7(0, P) has a rich structure which can be tested in the ARPES experiments. In observable quantities, 
the momentum dependence appears twice: via the matrix element A(P) and via the action 5(0) => 5(0, P). The 
analysis is simplified for the regime 2: the low polaron boundary Wi. Here the action dependence on O and P 
comes through the single variable O O + P 2 /2M\ where Mi ~ toAq/wq, is a heavy mass of the polaron center 
motion. This kinetic energy contribution can be neglected in compare to the matrix element dependence on P which 
confines A 2 <~ \^p\ 4 within the characteristic momenta distribution \^p\ 2 of the wave function \P(a;) of the self- 
trapped electronic state localized over the scale L ~ £ ; beyond P ~ ^q 1 = A /hvp, the function A(P) falls off 
exponentially. (At this scale, the recoil kinetic energy P 2 /2M ~ uj 2 /Aq is small in compare to the energy width 
e ~ (5") -1 / 2 ~ (u>q(W 1 — 0)) 1//2 ^> loq. Then the final integration over P affects only A(P) and gives a constant 
factor ~ l/£o-) Altogether we find for tunnelling just the law @ with O — ► U/2. 
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In the regime 1., near the free edge, the states are shallow e = 1 — E/Aq << 1 and extended L/^q ~ e x / 2 3> 1. 
The effective mass M for the center of motion becomes light, energy dependent 

e 2 e 5 / 2 

M~ -j- r , 

w 2 L w 2 

but the characteristic energy scale of the form factor ~ M _1 L -2 ~ oj 2 e~ 3 ^ 2 is still small in comparison with the 
characteristic energy width ~ J^ 3 of Q. So again we integrate separately the factor J A 2 (P)dP ~ J dP^j^ ~ L to 
obtain an additional prefactor e -1 / 2 for the tunnelling law JTJ with £1 — * {7/2. 

Recall that for the ARPES with independent variations of f2 and P, their interference may lead to rather unexpected 
and potentially observable phenomena (4, section III.D). One of them is the "quasi spectrum" : the intensity maximum 
over the line Ao — fi ~ P 1 / 2 within the PG Q < Ao (-4, section III.D, case Bl, Eq.49). Another effect is the emergence 
of instantons at high P within the domain of free electron region > fl > Ao leading to the enhanced intensity within 
the band < Q - A < cnst/P 3 (£, section III.D, case B3, Eq.51). 

III. TUNNELLING: THE DERIVATIONS. 

We shall follow the adiabatic method of earlier publications^^ assuming a smallness of collective frequencies ujq 
in compare with the electronic gap: ojq -C Ao. Now, electrons are moving in a slowly varying potential A(x,t), so 
that at any instance t their energies Ej(t) and wave functions ^j(x,t) are defined from a stationary Schrodinger 
equation H^(x,t) = E(t)^(x, E(t)) (Eq. (|A1|) below will give an example). The Hamiltonian H = H(x,A(x,t)) 
depends on the instantaneous configuration A(x,t) so that E(t) and ^f(x,E(t)) depend on time only parametrically. 
Exponentially small probabilities which we are studying here are determined by steepest descent paths in the joint 
space [A(x), t\ of configurations and the time, that is by a proximity of the saddle point of the action S. It is commonly 
believed, in analogy with the usual WKB, that the saddle point, the extremum of S over A and t, lie at the imaginary 
axis of t so that, as usual, we shall assume t => it and correspondingly S iS since now on. 

Consider the system of two weakly coupled chains j — a,b which are put at the electric potential difference U. The 
system is described by the total action 

S ab = S a + S b + t 1 _ [ dxdt{i>i{x, t)%(x, t) + fy(x, t)#.(x, *)), (5) 



where Sj = S[Aj(x, t)] are single chain actions and the term ~ t± describes the interchain hybridization of electronic 

sates. &j(x) are operators of electronic states. 

The average transverse current is given by the functional integral 

= J D[A 3 (x, t)] it ± (& a (x,t)y b (x, t) - ft(:c, t)* a {x, t)) exp[-S ab ] 

JdxDiA^x^t^eM-Sab] ' 1 ' 

A. One electron tunnelling. 

We consider first the processes originated by the transfer of one electron between the chains. They appear already 
in the first order of expansion of the exponent in JSJ in powers of t± , which contribution to the current JfJJ can be 
written as 

J = Z^ 2 tl J D [A, (x, t)] J d(x -y) J d( n - t 2 ) (7) 

[*l(x, n)* fe (x, n)V a (y, T 2 )V* b (y, r 2 ) exp(-5(n - r 2 , Aj{x, t))), 

where the normalizing factor Zq 2 is the denominator in © taken at t j_ = 0. Here the time dependent action S(ti — r 2 ) 
describes (in imaginary time) the process of transferring one particle from the doubly occupied level E a < of the 
chains a to the unoccupied level Eb > of the chain b at the time t\ and the inverse process at the time r 2 . We have 

/ +/ \dt[L a (0) + L b (0)}+ dt[(L a (-l) + L b (l)-U] 



DC 



dt[L a (0) + L b (0)} + / dt[(E b +E a - U], (8) 
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where Lj{v) = L([Aj],i/) are Lagrangians of the j-th chain with the number of electrons changed by v. They are 
given as a sum of the kinetic term and the potential V v : 

Lj(y) = [ dx^¥- +V v [A(x,t)] ; V v = Vo + \v\E. (9) 

Here the potential term V u contains the energy of deformations and the sum over electron energies in filled states a 
which include both the vacuum states and the split off ones: 

V v [A(x,t)]= fdx^+ ]T E a [A(x, t)] — Was (10) 

J 9 E a <E F 

(here g is the coupling constant). V v is counted with respect to the GS energy Was so that in the non perturbed 
A = Ao state V u — \v\Aq (the particle, electron for v > or hole for v < 0, added instantaneously to the non 
deformed GS is placed at the lowest allowed energy, the gap rim Ao). 
The exact extremal (saddle point) trajectory is defined by equations 



SS/6Aj(x, t) = , dS/dn = dS/dr 2 = 0. (11) 

Actually the explicit calculation of the action requires for approximations. We shall follow a way^ of the zero 
dimensional reduction which reduces the whole manyfold of functions Aj(x,t) to a particular class 

Aj(x,t) A E (x-Xj(t),Ej(t)), S[Aj(x,t)} => S[E,(t),X 3 (t)} (12) 

of a given function A E oi x (relative to a time dependent center of mass coordinate Xj(t)). A(x) is parameterized 
by a conveniently chosen (see^ for examples) parameter for which a universal and economic choice is the eigenvalue 
Ej(t). The requirement for the manyfold A(x,E) is that it supports a pair of eigenvalues ±E split off inside the 
gap (— Ao,Ao) which span the whole necessary interval. The last simplification is to assume, in the spirit of all 
approaches of optimal fluctuations^*^, that the potential supports one and only one pair of localized eigenstates 
^(x, ±E). Explicit formulas for the Peierls case are given in the Appendix El 

Recall that for the OA problem we deal with one chain characterized by one pair of functions E(t) and X(t). But 
for the interchain tunnelling, the functions Ej(i) at chains j = a,b are not obliged to be identical and also the wells 
may be centered around different points Xj(t). Within such a parametrization the variational equation in i|ll|) yields 
the equation of motion for E(t) 

m){^fj -V^-^U + H^ = 0, f(E) = -^Jdx(^^iy , (13) 

where H V j = cnst are the Hamiltonians which must be constants within each interval of integration in (jSJ . Apparently, 
at the outer intervals (t < t\),(t2 < t) Hoj = to provide the return to the GS with Vq = at t — > ±oo. At the inner 
interval (t\ < t < T2) H\j = Ej(ri) + U = Ejfa) + U to preserve the continuity of velocities Ej at t = n,2- Since 
the values Ej (ti^) are determined uniquely by the equation of motion at the outer intervals, then Ej(ri^) coincide 
for both j = a, b, hence H a = Hi, and the functions Ej(t) become identical at any time E a (t) = Ef,(t) = E(t). (Still, 
the shapes are allowed to be shifted by different centers Xj(t): A a (x — X a ,t) = Ai,(x — Xb,t)). Finally the extremal 
conditions Ijlll) with respect to impact times t; in (JSJ yield 

E a (t) + E b {t)\ t=Ti2 = U hence E a { Tl<2 ) = E b {r ia ) = U/2. 

The action is finite S < 00, hence the transition probability is not zero, only for a closed trajectory, that is at 
presence of a turning point (as examples, sec figures I5IKI7IXI in the Appendix EJ. There must be a minimal value of 
E = E m where E = hence V\{E m ) = U/2 and E m < U/2. The last condition requires for minV\{E) =W\ < U/2 
that is for U > 2W\ which determines the threshold voltage at twice the polaron energy. 

We arrive at the effective one chain problem with the doubled effective action. The extremal tunnelling action 
is Stun = 2Si which is twice the exponent appearing in the spectral density I with limiting laws (llll't . The full 
expression is 

rU/2 /.Ao 

Stun(U) = 8 / dEy/fiVj. - U/2) + 8 / dE^/fV ; Vx{E m ) = U. (14) 

JE m JU/2 
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We obtain a final expression for the current after integration over Aj(x, t) around the extremal taking into account 
the zero modes related with translations of the instanton centers positions Xj (t) . (Details of calculations are given in 
the Appendix lB|) 



J(U)KtlMuLu a ^j J ^e-" 2i2 / 4 |vl/ f) ((7/2)* p (-C//2)| 2 cxp[-25 / ([//2)], 



(15) 



where ty p is the Fourier transforms of the wave functions ^(x), the time T is defined as T = JjjJ^ dE/E. The mean 
fluctuational displacement / of the center of mass between the impact moments is given as 



i 1 -- 

where M{E) is the translational mass: 

M(E) 



dt 



M(E(t)) 



dx 



u ' 2 dE 



VfW) 



M(E) y/ViW-Uft' 



dA(x,E) 
dx 



Mu = M{U/2). 



(16) 



Note that the prefactor in Eq. 1)150. which is the matrix element between orthogonal states ^f(E) and ^(—E), is 
always nonzero due to the integration over zero modes Xj (t) (in contrast to results for the rigid lattice where it obeys 
the selection rules); see more in the Appendix iBl 

Comparing with the PES intensity 7(0) calculated in^ we see that, up to pre-exponential factors, the tunnelling 
current is proportional to the square of the PES intensity I: J cx t\l 2 (£l — U/2). E.g. near the threshold U — 2W\ 
we can write 



J 



Aoia 



Ao 

L0() 



3/2 



exp 



exp 



(U - 2Wi) 

C 2 log 

guj 



2C 3 A 



{U - 2Wi) 



(17) 



The coefficients Q ~ 1 can be found numerically from (|14l) as (for the Peierls model) C\ = 0.4, C 2 = 2.9, C3 = 0.1. 
(These values differ from the corresponding ones in£ because of different normalizations of frequency ujq in compare 

tO OJph)- 



B. Bi-electronic tunnelling. 



It is known that the joint self-trapping of two electrons allows to further gain the energy resulting in stable states 
different from independent polarons. In general nondegenerate systems this is the bipolaron, confined within the 
length scale twice smaller than that of the polaron, the energy gain of the bound state SE = Ao — E is four times 
that of the polaron and the total energy gain of the bipolaron 8W2 = 2Ao — W2 is also four times that of two 
polarons. (Certainly these results neglect the energy loss due to the Coulomb repulsion which may become critical 
for the stability of a shallow bipolaron.) The same time, the total energy of one bipolaron Wi — 2Ao — 6W2 is larger 
than the energy of one polaron W\ = Ao — 5W\ and even than the free electron energy Ao. This is why bipolarons 
cannot be seen as thermal excitations while they are favored in case of doping. The information on their existence 
comes from the ground state of doped systems where bipolarons are recognized by their spinless character and special 
optical features (see 14 for experimental examples on conducting polymers and 15 for relevant theoretical models). An 
important advantage of tunnelling experiments is a possibility to see bipolarons directly, at voltages U below the 
two-polaron threshold 2W\ that is within the true single particle gap. This possibility comes from the fact that, for 
bipolarons as particles with the double charge 2e, the voltage gain by transferring from one chain to another is 2U, 
hence the threshold will be at U — W 2 < 2W\. The probability of the bi-electron tunnelling is small as it appears 
only in the higher order ~ t\ in interchain coupling. But it can be seen as extending below the one-electron threshold 
where no other excitations can contribute to the tunnelling current. 

The bi-electronic contribution to the current can be written, by expanding (J5J and JBJ, as 

J 2 = Z^ 2 ti f D[A a ]D[A b ] [fld yi dnexp(-S{n)) 
J J i=l 

[*S(x,T)* i (l,T)**(ltt ) r 1 )* i (tfi,Tl)*J(yi,75)*a(lfa S 7a)*;(»3 ) 7s)*B(|ft,7a) -{2 «-> 3}], 
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which generalizes expressions (7) and (8) for the one electron tunnelling. Here 
S({n}) = jy +y }dt[L a (0) + L 6 (0)] + 



T3 

T3 ~| I-T2 



dt[(L a (l) + L b (l) - U] + J dt[L a (2) + L b (2) - 2U}. (18) 
Within our model (|9I10() the potentials V are additive in energy E, then the action can be simplified as 



S({n})=2 dtL a (0,t) + 

J —OQ 

dt(E a (t) + E b (t) - U) + 2 f 2 dt(E a (t) + E b {t) - U). (19) 



Tl 



The extremal solution is defined, as above, by equations of the type (|11|) but with four impact times Tj instead 
of two. (Actually, in view of the time reversion symmetry, the number of boundary conditions is twice smaller.) 
A similar analysis of the extremal solution shows that optimal fluctuations A.j(x,t) are identical in shape, up to 
shifts of their centra: Xj(t), A.j(x,t) = A(x — Xj,t). Hence the energies are identical E a (t) = E b (t), and also the 
resonance conditions 2E(ji) = U take place at the impact moments Tj. Moreover, the simple hierarchy of our model 
V2 — V\ = Vi — Vo = E — U/2 shows that all branches V V {E) cross at the same point E = U/2 (see figures 15^6171 
below). Then the evolution E(t) switches directly from the branch v = to the branch v = 2 and back, without 
following the intermediate branch v = 1. It means that the intervals (r, ti) and (t2,T3) of one-electron transfers 
v = 1 are confined to zero: t = t%, t% = T3. In other words, only processes of simultaneous tunnelling of pairs of 
particles are left. Notice that this picture changes in more general models, particularly taking into account important 
Coulomb interactions. They add, to the energy branch of a shallow bipolaron, the energy 6V 2 ~ (e 2 /e±L) \n(L/a±) 
where e± is the dielectric susceptibility of the media in the interchain direction, L = L(E) is the localization length 
of ^(x, E), such that E ~ l/(mL 2 ). Now the intermediate intervals (r = ti),(t 2 = r 3 ) appear where the evolution 
follows the v — 1 branches, see figure |H1 With increasing Coulomb interactions this single particle interval becomes 
more pronounced and the bipolaronic threshold is shifted towards the one of two independent polarons. 

In any case, the extremum solution for the action i|18f) is achieved on the instanton trajectory given be the equation 

f(E)E 2 = Vu(E) = min{foCE),Vi(£) - U/2,V 2 (E) - U}. 

The extremal action is 

/•Ao 

S 2 {U) = 8 / dE^/f(E)Vu{E) ■ Vu{E m ) = V 2 (E m ) - U = 0. (20) 

J E m 

This action is finite if the turning point E m does exist, that is if U > minV2 = W%. 

Notice that, neglecting Coulomb interactions, the energy V v is determined only by the total number v = v e + v b = 
u{E) + (2 — u(—E)) of electrons and holes. Then the energy of the bipolaron (both v{E) and v(—E) are cither empty 
or doubly occupied) and the energy of the exciton (both v(E) = 1 and v{—E) = 1 are singly occupied) are the same. 
Then the trajectory of the bi-electronic tunnelling becomes the same as the one for the case of optical absorption 3 , 
only the action is doubled S 2 (U) = 2Soa(Q — U). Up to the pre-exponential factor we have 

j 2 Kti[i OA (n = u)] 2 , (21) 

where Ioa(Q) is the optical absorption probability for one chain. 

For common systems with a nondegenerate ground state, the dependence S 2 (U) resembles qualitatively the law !17l for 
the one-electron contribution, with a similar behavior near the threshold U — 2W\ — > U — W 2 . The situation changes 
for a doubly degenerate ground state where the bipolaron dissolves into a diverging pair of solitons (dimerization 
kinks). Thus for the Peierls model the evaluation of l|2U|) gives, similar to the OA law of 3 -, near the two particle 
threshold 



J 2 ~ *i exp - max S 2 H J6A Q (U - 2W S ) ; max5 2 = C4A0/5W0 (22) 

V 9^o J 

with C4 = 3.77. The overall dependence for the log J 2 {U) ~ —S 2 (U) is shown at the figure [21 Here we see explicitly 
that in the order ~ t A ± the threshold voltage U — 2W S is smaller than U = 2W P obtained in the order ~ t\. Therefore 
this is the main contribution to the current in the region 2W S < U < 2W p . Figure |3 shows that the dependence 
of J 2 (U) near the low U onset is much sharper than that of J(U) at the figure ^ near the polaronic onset which 
corresponds to the higher singularity in the limiting formula 1221) in compare to l|17|) . 
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IV. DISCUSSION AND CONCLUSIONS. 

In quasi ID systems with a gapful electronic spectrum, the interchain tunnelling (as well as PES or OA) can be 
used to test virtual electronic states within the pseudogap. Due to the interaction of electrons with a low frequency 
mode, phonons in our examples, the tunnelling is allowed in the subgap region U < E® which forms the pseudogap. 
The one electron processes lead to universal results similar both for systems with the build-in gap and for those where 
the gap is due to the spontaneous breaking of a discrete symmetry. The PG is entered with the law Q and continues 
down to the threshold U\ = 2W\, approached with the law J5J). This threshold corresponds to the interchain transfer 
of fully dressed particles: polarons with the energies W\. But in tunnelling the PG is stretched even further down 
thanks to processes of a simultaneous tunnelling of two electrons. It terminates at the lower threshold U2 = W2 or 
U2 = 2W S , U2 <U\. Here W2 is the energy of the bipolaron - a bound state of two electrons selftrapped together. In 
degenerate systems the bipolaron dissolves into unbound solitons, hence the threshold at 2W S with a more pronounced 
dependence of the tunnelling rate (1221 as well of the OA. Numerical results are presented at figures HTH 

There is an important difference between subgap processes and the usual overgap transitions at U, fl > E® of free 
electrons in a rigid system. It comes, beyond intensities, from different character of matrix elements. Actually within 
the PG region there are no particular selection rules since the wave functions of virtual electronic states split off 
within the gap are localized having a broad distribution of momenta. Then the PG absorption is allowed independent 
on the interchain ordering. Contrarily, the regular tunnelling across the free gap shows an expected DOS singularity 
~ {E—E®)~ x / 2 for the out of phase interchain order while for the in-phase order the threshold is smooth ~ (E— Eg) 1 / 2 . 
This difference may be important to choose an experimental system adequate for studies of PGs. The smearing of 
the free edge singularity is a natural criterium for existence of the PG below ibi&. But the total absence of this strong 
feature in systems with forbidden overgap transitions can allow for a better resolution of the whole PG region, down 
to the absolute threshold. Probably a very smooth manifestation of gaps in usual tunnelling experiments on CDWs 16 , 
while the gaps show up clearly through activation laws, is related to this smooth crossover from the overgap to the 
subgap region. (Notice that the existing experiments refer mostly to ICDWs which, with their continuous degeneracy 
of the GS, must be studied specially which is beyond the scope of this article.) 

Finally we shall discuss relations with other theoretical approaches. Most theories of tunnelling, see£, keep the 
following assumptions: i. They refer to the overgap region where interactions or fluctuations are not important and 
usually are not taken into account, ii. They refer to the incoherent tunnelling, local in space, which is a usual 
circumstance of traditional experiments. The PG in tunnelling was considered by Monz et al in 7 in the framework 
of the approachi 7 -. This method became popular recently in theories of the PG thanks to its easy implementation: 
it is sufficient to average results for a rigid system over a certain distribution of the gap values. Apparently this is 
the way to describe an average over a set of measurements performed on similar systems with various values of the 
gap, e.g. manipulating with the temperature, the pressure or a composition. But actually, as we could see above, the 
PG is formed by fluctuations localized both in space and time, the instantons, with localization parameters depend 
on the energy deficit being tested. There is an intermediate approach applied 1 ^ to a complex of the PG phenomenon 
from optics to conductivity and susceptibility. It treats fluctuations as an instantaneous disorder due to quantum 
zero point fluctuations of the gap. Indeed, this picture can be well applied, as it was done already in2i, but only to 
dynamical processes and only in the upper PG region, just below the free gap E°, which leads to the law QJ. But 
deeper within the PG, the fluctuations are not instantaneous: they require for an increasingly longer time and become 
self-consistent with the measured electronic state leading to another law and to appearance of the lower threshold. 
Generalizations and deeper analysis of the model of the instantaneous disorder lead to interesting theoretical studies^, 
but their applicability is very limited unless the variable time scale is realized as we have demonstrated in this and 
preceding articles. 

Our approach can be compared to the work«2£ on the fluctuational creation of pairs of phase solitons in a ID com- 
mensurate CDW under the longitudinal electric field. But in our case me deal, in effect, with the interchain tunnelling 
of pairs of solitons under the transverse field; also the solitons have a more complex character of a multielectronic 
origin. 

In conclusion, the presented and earlies^^ studies recall for the necessity of realizing the variable time scale of 
subgap processes both in theory and in diverse interpretations of different groups of experiments (dynamic, kinetic, 
thermodynamic) which address excitations with very different life times. 
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APPENDIX A: DETAILS ON SELF-TRAPPING BRANCHES. 



We consider the system of weakly coupled dimcriscd chains. Each chain is described by a usual electron- phonon 
Hamiltonian (Peierls, SSH models). Electron levels E and wave functions \& = &(x, E) are determined by equations 

[-v F id x a 3 + A(a)o-i]* = EV, (Al) 

where 01,3 are the Pauli matrices, = (^+, \&_), ^±(2:) are the components of electron wave functions near Fermi 
points ±Pf j and the real function A(x) is the amplitude of the alternating dimerization potential. The ground state 
of each chain is the Peierls dielectric with the gap 2Ao. The electron spectrum has the form E% = v 2 F p 2 + Aq (in the 
following we shall put the Fermi velocity vp = 1, A = 1 and, as everywhere, the Plank constant H = 1). The excited 
states are solitons (kinks), polarons and bi-solitons (kink-antikink pairs) which are characterized by electron levels 
localized deeply within the gap (see the review^). The one parametric family of configurations A(x, E) supporting 
the single split-off pair of levels ±E can be written as 

2(1 — E 2 ) 

As, 2-0=1 K 7 1 77^ (A2) 

l + £:coshf2x(l-£; 2 ) 1/2 J 

evolving from a shallow potential well at E « 1 through the stationary configuration for a polaron v = 1 to the pair 
of diverging kinks at E — ► as shown at the figure [3] The potentials V v (for the E level filling v) as functions of E 
are given as 

VJE) = vE+-\?l-E 2 - -Ecos^E. (A3) 

TT TT 



The translational mass can be found as 
M(E) = 



- tan 3 [ cosh" 1 — ) - E 2 cosh" 1 — + E 2 y/l - E 2 
3 V E E 



(A4) 



Consider the matrix element between levels ±25 in the Peierls state. The wave function has two components (u, w) 
according to \& = ^(x, E) — u cosp^x + w smppx. Explicit expressions for split-off states are u,w ~ y/l — A 2 ± d x A. 
The equation for the bound eigenstate (|A1|) shows the following symmetry: w(x, E) = u(—x, E), w(x, —E) = u(x, E), 
u(x,—E) = —u(—x,E). Then 9(x, E) — ((u(x),u(—x)), ^>(x,—E) = (—u(—x),u(x)), with u — u(x,E), which 
demonstrates explicitly the orthogonality of ^(x,E) and *&(x, — E). The matrix element in Eq. (|15J) becomes A 2 ~ 

\^p(E)'9p(— E)\ 2 = \— u p u-p + u p Up\ 2 . At p = 0, A = Ao = 0, hence for identical chains the transition at the free gap 
fl = 2Ao is forbidden which removes the singularity at the gap threshold in a rigid system. But the true threshold 
at 2Wi for the subgap absorption or tunnelling are not subjected to this selection rule since the wave functions of 
localized states associated to the optimal fluctuation are distributed over the momentum region p ~ CcT 1 - 

Figure 13 shows exact shapes A(x, E) of the equilibrium polaron (upper thick line) and of a well formed (E = 0.01) 
pair of solitons (lower thick line). Thin lines show exact shapes of optimal fluctuations necessary to create these states 
by tunnelling. Notice the much less pronounced shapes for optimal fluctuations in compare to the final states which 
facilitates the tunnelling. 

Figure 0] plots the total energy V\ (E) of the single particle branch as a function of the associated energy of 
the bound state. E = 1, Vi(l) = 1 corresponds to the particle added to the unperturbed ground state, at the 
bottom of the continuous spectrum. E = is the mid-gap state reached for the limit of two divergent solitons when 
the total energy approaches the maximal value Vi(0) = 2W S — A/tt w 1.27. In between, at E\ = 2" 1 / 2 sw 0.7, 
V\{2~ 1 ^ 2 ) — W\ — 2 3 / 2 /tt ps 0.9, the minimum corresponds to the stationary polaronic state. The short thin vertical 
line between plots E and V\(E) points to the configuration (upper thin curve at the figure [3} of the fluctuation 
necessary for tunnelling to the polaron (the minimum of V\ (E) , upper thick curve at the figure 2J) . 

Next three figures plot the total energies V v — vU for different branches as a function of the energy E of the associated 
bound state (all in units of Ao). Branches are distinguished by their ordering at E = 1. Figure [5] corresponds to the 
potential U = 1.2 which is below the bi-electronic threshold; no branch is crossing V — axis, hence no final action 
is allowed and the current is zero. 

Figure corresponds to the potential U = 1.4 which is between the bi-electronic threshold 2W S = A/tt w 1.3 and 
the polaronic one 2W\ — 1.8; the bi-electronic branch crosses the axis V — at the point E m , the action is finite, 
hence a nonzero tunnelling of two electrons is allowed. 

Figure [7] corresponds to the potential U — 1.8, above the bi-electronic threshold 2W S — A/tt w 1.3, exactly at the 
polaronic one 2W\ = 1.8. Now two parallel processes of one- and two- electron tunnelling are allowed. 
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FIG. 3: Exact shapes A(x,E) of the equilibrium polaron, E — 2 -1 / 2 (upper thick line) and of a nearly formed (E = 0.01) 
pare of solitons (lower thick line). Thin lines show exact shapes of optimal fluctuations necessary to create these states by 
tunnelling. 
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FIG. 4: Total energy of the single particle branch Vi (E) as a function of the energy E of the associated bound state. 
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FIG. 5: Total energies V v — uU for different tunnelling branches V v — vU as a function of the energy E of the associated bound 
state. This figure corresponds to U = 1.2 which is below the bi-electronic threshold. 
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FIG. 6: Tunnelling branches V u — vJJ as a function of the energy E of the associated bound state. This figure corresponds to 
U = 1.4 which is between the thresholds for tunnelling of bipolarons and polarons. 




FIG. 7: Tunnelling branches V v — uU as a function of the energy E of the associated bound state (over a selected interval). 
This figure corresponds to U = 1.8 which is just at the polaronic threshold, above the bipolaronic one. 

The figure |H| corresponds to the potential U = 1.6 between the bi-electronic threshold 2W S — 4/tt ps 1.3, and the 
polaronic one 2Wi — 1.8. Contrary to the figure the Coulomb interaction is taken into account which lifts the 
degeneracy of the earlier crossing point of three branches. The one electron term v — 1 does not cross V = axis 
yet, but it passes below two other terms in a vicinity of their crossing. Now the optimal bi-electronic tunnelling takes 
place via a sequence of two single electronic processes confined in time. 

APPENDIX B: DERIVATION OF PREFACTOR 

We need to perform the integration over A.j(x, t) around the extremal taking into account the zero modes related 
to translations of positions Xj(t) of the instanton centers. The path integration over the gapless mode X(t) is 
important, particularly for the matrix element: the overlap of wave functions evolves following X a {t) — Xb(t) while 
their localization follows the evolution of E(t). We shall work within the zero dimensional reduction of Eq. <|12[) . 

We expand the field Aj(a;, t) in the vicinity of the instanton solution as 

A j (x,t)=A (x-X j (t),E j (t) + S(x-X j (t),t)). (Bl) 

Following-, we rewrite (JSJ) as 

J K *i II / d ( x - yM T i - T2 ) / D[Ej\D[Xj]J Xj Je 3 exp(-5) 

j=a,b 
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FIG. 8: Tunnelling branches V u — vJJ as a function of the energy E of the associated bound state. This figure corresponds to 
U — 1.6 which is between the thresholds for tunnelling of bipolarons and polarons. Contrarily to the figure^ the Coulomb 
interaction is taken into account which lifts the crossing degeneracy. 



- X a ( n ), E{ Tl ))^ a {y - X a (T 2 ),E(r 2 ))* b (x - X b ( n ), -E(n))*t(y - X b (r 2 ), -E(r 2 )), (B2) 

where Jx =oc Iln=i y/M(E(t n )), Je, oc n n =i \/ f(-®j(t n )) are the Jacobians of the transformation lfT2*)l . (N — > oo is 
the number of points for the intermediate discretization of the time axis.) We integrate over the zero mode X(t) and 
take into account fluctuations of the instanton shape due to variations of the parameter Eo(t). 
The action in (|12fl has the form 

S[X, E] = ]T dt (M(E 3 (t))X 3 2 /2 + m)E]/2 + V^E-)) 

j=a,b 

with Vu (E) from (|2(J[I . The integration over DXi (t) is carried out exactly after the transformation MX 2 = Z 2 using 
the known expression 



d 2 S r 



D[x] exp ( - / dt{— + V(x)) ) ~ C M-Sci)]lj—^, (B3) 



where x\ = x(ti), x 2 = x(t 2 ). Next, we perform in (IB2f) the remnant integrations over coordinates at the impact 
moments: X x = X a (n), X 2 = X a (r 2 ), Y x = X b (n), Y 2 = X b (r 2 ): 



Jocti / dxdT X dX Y dX 2 e V"-**^ - Xx, E{t x ))^ a {y - X 2 , £(t a ))2^_I 

dY x dY 2 e ^%{x - Y u -E{ti))* b {y - Y 2 , -£(*.,)) X_i exp[-S E (E)]. (B4) 

h 

Here Mj = Mj(ri 2)1 and the same for l i; are functions of energies in these points which finally become E^ = Ei(ri 2 ) — 
U/2. Using Fourier transforms, we rewrite the product of wave functions as 

Integration over dxdy gives S(p — p), 5(q — q) and integration over dX 2 , dY 2 gives L8{p — q). After integration over 
d(Xi - X 2 )d(Yi - Y 2 ) we arrive at the result (JT5). The factor ^JdT/dU = l/yJd 2 S/dT 2 in (JT5J after integration over 
D[Ei(t)] which was performed using again the equation l|B3(l . 
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